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' A microscopic model which has proven useful in describing amphiphilic aggregates as inhomo- 

H ■ geneities of a fluid is extended here to study the case of a two component surfactant mixture. We 

_ ■ have chosen an effective interaction between the amphiphiles that mimics the mixture of cationic- 

anionic surfactants. In agreement with experiments, and other theoretical approaches, we find 

■ regions where spherical vesicles are stable, with a well defined radius. The experimental dependence 
of the radius on the chemical environment is also qualitatively reflected in our model. 
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, , I. INTRODUCTION 

The structure and thermodynamical properties of fluid membranes composed of amphiphilic molecules are of prime 
scientific and technological importance. In particular the formation of vesicles has received very close attention 
due to their pharmacological applications (microencapsulation) and their use as models for biological membranes. 
G Usually, large vesicles are formed by intense mechanical shear (sonication) , but they are a metastable phase which 
disappears after a time of hours or days. Only recently has been found spontaneous formation of long-lived vesicles, 
suggesting that the vesicles could be an equilibrium state of the system (lj-|6| . We considerer here binary surfactant 
mixtures, which are the main body of the experimental systems with this behavior although there are also a 

few observations for systems with only one surfactant component |||| . Note that the thermodynamic status of these 
systems is still an open question. The spontaneous vesicles may be truly at equilibrium or they may be a long-lived 
transient state to a lamellar phase j?]]. It has also been argued that the existence of vesicles, and their size distribution, 
may result from the edge energy and the dynamics of growth . Else, as an equilibrium property, that they arise from 
the entropy associated with the position and shape of the vesicles || , rather than to the existence of a preferential 
curvature for an isolated vesicle. 

■ We may distinguish two main kinds of two-component systems in which spontaneous formation of vesicles appear. 
| In the first, the two amphiphiles have different self-assembly behavior. One of them (a phospholipid molecule) is 

quasi insoluble in water and tends to form lamellar phases, while the other (a surfactant molecule) is soluble in 
water up to a critical surfactant concentration and then tends to form micelles. One example of these systems is the 
mixture of octylglucoside and phosphatidylcholine [Q. The second possibility is that both surfactants have similar 
self-assembly properties and its vesicular behavior is due to the nature of the interaction between the two species. 
The main examples of these systems are mixtures of two ionic surfactants, with oppositely charged head groups. The 
mixture of the two types, in the solution, may form spontaneous vesicles, which would not appear in pure solutions of 
either of the surfactants. Many experimental observations have been reported on these catanionic systems ^-H]. In 
this work, we will only study this later kind of vesicular system, and search for the possible existence of spontaneous 
curvature, and typical vesicle sizes, based on the thermodynamic stability of single vesicles. We will neglect the effects 
of polydispersity and shape entropy. We will use the words vesicle and membrane indistinctly to refer to bilayer 
aggregates; in each case, we will make clear whether we are referring to a closed vesicle or a planar membrane. 

Recently Safran et al |lO|JTl| ] have treated this problem under the hypothesis that the existence of spontaneous 
curvature for the vesicles is associated with the breaking of the composition symmetry between the two molecular 
layers in the membrane. In a semiempirical Landau theory, the free energy of a spherical membrane of radius R is 
written in terms of its curvature, C = 1/R, and the relative concentrations of each species, A and B, in each molecular 
layer, xf = 1 — xf , with i = 1,2, for the two layers of the bilayer membrane. These relative concentrations are 
described by an overall relative concentration 

i> = + x£-xf- ) = i(a£ +x£- 1), 

which takes values ±1 for pure membranes (of A or B molecules) and is zero for a membrane with symmetric 
concentrations of A and B. The second order parameter is a transverse relative concentration, 

<!>=\{x$-xt) = \{xi-^)> 
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which measures the asymmetry between the two layers in the membrane, being zero for symmetric membranes. Since 
all the x\ have to be between and 1, (f> is constrained to be between 1/2— and — 1/2+ \tp\. For planar membranes, 
the free energy should be symmetric in <f), but there may be spontaneous symmetry breaking when the membrane goes 
through a phase transition which preferentially segregates the two species on different sides of the bilayer membrane. 
For vesicles with C / there is a difference between the inner and the outer layer, which breaks the symmetry of 
the free energy with respect to cj>. The presence of a coupling term proportional to C<j) favors the formation of curved 
vesicles with asymmetric composition C oc <f>. The resulting phase diagram exhibits a line of second-order transitions 
between the lamellar and vesicle phases, with the highest critical temperature corresponding to the system with a 
symmetric global composition i/j = 0. At lower temperature, there is a tricritical point where the transition becomes 
first order. Although these semiempirical theories give a very interesting and clarifying picture of these systems, they 
are parameter dependent and have a simplified functional dependence of the free energy. To check the validity of the 
later and to relate the parameters to molecular interactions requires microscopic theories. Safran et al Jl(]|ll| discussed 
the use of their model to explain the energetic stabilization of two-component vesicles, for mixtures of amphiphilcs 
which differ only in their polar head groups. The curvature is associated with the different concentrations of the two 
components in the inner and outer monolayer of the membrane. Also, they assume that this different concentration 
between layers is promoted by a interaction which gives phase separation in the same layer. As a consequence, they 
choose a repulsive interaction between the two different components. 



We have recently developed a simple microscopic model [12 13] for the study of aggregates of monocomponcnt 



amphiphilic molecules in water, which makes possible a consistent thermodynamic treatment of the system by taking 
account of the coexistence with the dilute solution of amphiphiles. The model considers only the positional and 
oricntational distribution of the amphiphiles, ignoring internal degrees of freedom. The effective interactions between 
the amphiphiles is described by a potential which includes the interaction with the water (so that water molecules 
need not be explicitly taken into account). This minimal model gives a qualitatively correct global phase diagram for 
aggregates with planar, spherical, and cylindrical symmetry. The elastic constants obtained are in good agreement 
with experimental results. Depending of the interaction potential, the model exhibits coexistence between a solution 
of free molecules, micelles, and planar membranes. This monocomponent model does not give spontaneous curvature 
for any interaction potential and therefore the spherical vesicles are always unstable when compared with the micelles 
or with the planar membranes. 

The goal of the present work is to extend the above model to the study of stability and phase behavior of two- 
component vesicles. As in the study of the elastic constants jt3|, the self consistent thermodynamic treatment, i.e. the 
thermodynamic equilibrium between the solution of amphiphiles and the vesicles, shows interesting restrictions on the 
possible values of the parameters used in semiempirical Landau theories. In order to exclude other possible sources 
of spontaneous curvature, we use a model in which the interactions are fully symmetric. Thus, the pure component 
membrane made with A molecules is exactly equivalent to the pure B membrane. In experimental systems with 
catanionic surfactants this symmetry is not perfect, i.e. the properties of the two types of pure component systems 
are not equal, but still they are much more similar than the properties of the surfactant mixture. 



II. DENSITY FUNCTIONAL MODEL 

In this section we extend our previous approach Jl^ ] to the case of two different components, which we label 
£ = A, B. The model is based on the description of the microscopic or mesoscopic molecular aggregates in water, as 
inhomogeneous density distributions, treated within the density functional formalism. We use an approximation to 
the exact Helmholtz free energy density functional, F[p^}, applied to a very simple molecular interaction model. As 
the amphiphilic molecules arc much larger that the water molecules, our previous approach described the interaction 
between them by a potential which represents an effective, water-mediated interaction ( and so, water molecules need 
not be explicitly taken into account in the model). The head part of an amphiphilic molecule attracts water much 
more strongly that the lipid tail. This fact produces an effective repulsion between the head of a molecule and any 
other molecule (regardless of the orientation of the latter because it displaces the water). We consider the amphiphiles 
as uniaxial molecules without internal degrees of freedom, other than the orientation of the head-tail direction, which 
is described by a unit vector u along the molecular axis. Let us consider a ^-molecule at ri, with orientation iii, and 
a ^'-molecule at r 2 with orientation U2. The effective pairwise potentials between them, with the proper features, 

may be written as (r 2 — ri, Ui, u 2 ) = <&q^ (r 2 i) + (r 2 i, Ui, u 2 ). The isotropic part $o only depends on the 
distance between the molecular centers, r 2 i = |r 2 i| = |r 2 — ri|, and the anisotropic part, $ a , may be written as a 
gen eral expansion in spherical harmonics. As we have argued in our previous papers on monocomponent systems 
[p2| , ^3| , we need only include the terms which do not couple the orientation of the two molecules. The general form is 
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(r 2 i, ui, u 2 ) = ( r 2i) x 

i=i 

[Pi (ui • f 2 l) + Pi (-U2-r 2 i)], (1) 

where f 2 i = r 2 i/r 2 i and Pi(cc) are the Legendre polynomials. This restriction of the model interactions has the 
advantage that minimization with respect to orientational degrees of freedom can be done analytically. 
The Helmholtz free energy density functional may be approximated by 



F[p A (r,ii),p B (r,u)}=F [p A (r),p B (r)] 
+k B T Y / dr i d UiP 4 ( r i) Q:C ( r i J Ui)log(47rQ! 5 (r 1 ,u 1 )) 



+ \ Y Y f dr 1 du 1 dr 2 du 2 p s (r 1 ,u 1 )pZ (r 2 ,u 2 )$| e (r 21 , u l5 u 2 ). (2) 

The first term is the density functional for the isotropic reference fluid; it depends on the density distributions, p*(r). 
The second term is the rotational entropy and the third one is the mean field contribution of the anisotropic interac- 
tions; both depend on the position and molecular orientation distributions, (fi(r, u), with a*(r, u) = p^ (r, u)//9* (r) . 

As stated in the introduction we are interested in the study of a mixture of surfactants with similar self-assembly 
properties, so that the formation of vesicles is due to the interaction between the two species. Within this scheme, 
it is logical to assume that the A- A interaction is equal to the B—B one, but different from the A—B interaction: 
$> AA — <& BB ^ <& AB . We have fixed the shape of the potentials following the experience gained in the monocomponent 

case. We consider that the isotropic potential energy (r) is essentially repulsive, since an attractive contribution 
would produce the unrealistic segregation of water-rich and amphiphilic-rich bulk fluid phases. As the simplest choice, 

we use a repulsive isotropic interaction consisting of a hard sphere (HS) potential, $^ (r). By symmetry, the sphere 
diameter, dhs, is taken to be identical for the two species. We choose this hard sphere diameter as the unit length and 
ksT = (i^ 1 as the unit of energy. For the anisotropic potential we have considered a single form, irrespective of the 
species involved, changing only their prefactors. The first coefficient function of the anisotropic potential given by eq 
1 was taken to be an empty-core Yukawa potential, 

te' 

$f (r) = — — exp(-\ (r - d hs )), 

for r > dhs, and null inside the hard core, r < dhs- The next coefficient function in eq 1 is taken to be proportional 

to the first, <j>| ? (r) = (r), and all the rest coefficient functions are neglected. Based in our previous work, 

we will consistently use Xdh s = 2 and q = 0.5. The latter choice is made so that the monocomponent planar 
membrane will be stable with respect to the micelles, in the range of couplings we consider. In contrast, when 
q < 0.3 the micellar structures are more stable. Thus, the model only has two unprefixed adimensional parameters: 
K = f3K AA /dhs = (3K BB /dhs, the measure of the amphiphilic character of the pure surfactants, and the ratio of 
like/dislike energies, given by e = Kab / Kaa- Note that within this model we could study the first kind of vesicle 
systems, in which the two surfactant have different self-assembly behavior, considering, e.g., q = for the surfactant 
component and q — 0.5 for the phospholipid molecules. This study is left for a later work. 

For a given temperature, T, and chemical potentials, p A ,p, B , we have to search for local minima of the grand 
potential energy, with respect to the distribution functions p^(r,u). Due to our simplified choice of anisotropic 
interactions, the minimization with respect to o^(r, u) may be carried out analytically. Thus, 



a«(r,u) = lexp I £ a«'u+ £ u T £«'u], 

\£'=A,B i'=A,B 



where the constants ensure the normalization condition of o^, i.e. J o^(r, u)du = 1. For each possible value of 

£ and £ , is a vector arising from the first coefficient function of the anisotropic potential and is a tensor 
associated to the second coefficient. Their values are given by: 



a«' (n) = -pj dr 2 /(r 2 )$f'(|r 12 |)?i 2 
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and 

3fi2 (g> fi2 



B«'(n) =-/3y dr 2 /(r 2 )$f'(|r 12 |): 



where I is the identity matrix and <S> is a direct product. For planar and spherical symmetries the expressions simplify 
and we only have to consider the modulus of (r) and the largest eigenvalue of S'* (r) ; these we will call (r) 
and (r), respectively. 

By substitution of the above into eq 2 we get a free energy density functional, which is already minimized with 
respect to the molecular orientations: 

F[p A (r),p B (r)} = mm{F[p A (r,u),p A (r,u)} } a = 

= AF hs [p(r)] + k B T I dripHri) log^fa) - 1) (3) 

£=A,B J 



I <V(r)log(Q«(r)), 

_ A D " 



£,=A,B 



where 



QHr) = \ 



1 \('=A,b , 



and p(r) = p A (r) + p B (r) is the total density. In our calculations the interaction part AFh s [p(r)] of the hard sphere 
free energy is approximated by a well tested non-local density functional [ff3 . 

In absence of external potentials, this density functional always has a stationary grand-potential energy for a 

homogeneous density distribution p^(r) = Pq, for which (r) and (r) vanish. This state represents a mixture 
of surfactant molecules with random orientations, which in our mean field description is equivalent to a hard spheres 
mixture. The densities of this homogeneous solution are used to control the chemical potentials, pfi, in our search for 
inhomogeneous density distributions. In this case the expression for the chemical potentials is 

= log(4) + (3p hs (p ), 

where Phs(po) is the excess (over the ideal) chemical potential of a hard sphere fluid at a total density pq. 

For low values of pfi (and p^), the homogeneous phase represents a dilute solution of amphiphilic molecules in water 
and we may find other local minima of the grand potential energy, 

n[p A (r),p B (r)} = F[p A (r),p B (r)} - p A J dv p A (r) - p B J dr p B (r), (4) 

for non-uniform density distributions, which represent the different types of molecular aggregates. 

III. RESULTS 

We are interested in the vesicle-planar membrane phase diagram. In this phase diagram we must require that 
the aggregates are in chemical equilibrium with the dilute surrounding solution of free molecules (which we call "the 
bulM'), i.e. that they be local minima of eq 4. Beyond of this condition of local stability, we have to seek a global 
stability, which means that the excess grand-potential energy of the aggregate, with respect to the bulk solution, has 
to be exactly zero. This is a strong requirement, which comes from the consistent thermodynamic treatment of the 
model, and it sets important restrictions on the possible structures of the aggregates. 

Although the minimization of eq 4 can be carried out exactly, we have calculated the coexistence line within a 
parameterized family of variational density profiles, which improves the computational efficiency without important 



changes in the results 12 1. With this approach, it is easy to impose a well-defined mean radius to the structure, whose 
inverse is the curvature, C. In order to search for vesicles, we impose spherical symmetry and minimize the free energy 
of our system with respect to the density profiles. As vesicles are characterized by two peaks corresponding to the two 
monolayers, a useful parameterization of each of the two density profiles is achieved by two Gaussians superimposed 
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on an uniform background; hence we will need four of them. The minimization is made with respect to the twelve 
parameters that characterize the Gaussians: widths oq (four parameters), surface densities of the monolayers rji (two 
parameters), relative compositions x\ (two parameters), radius (two parameters), and distances between layers S 
(two parameters). In practice, the radius and the distance between layers are almost always equal for both species. 
Thus, our density distributions are given by: 



1/2 




cxp[-a«(i?-(i?«+d«)) 2 



(5) 



Here, the x\ have the same meaning as in ref 10 and 11, and from them we get directly the variables <j) and ip defined 
above. 

In Figure 1 we show a typical profile obtained with our model for a spherical vesicle. As can be seen, the profile 
exhibits two layers of amphiphiles with different composition. In this case the concentration of the minority component 
B is higher in the outer layer than the inner one and the majority component has an opposite behavior. The 
orientational profiles are very similar to those obtained in the monocomponent case, being saturated and of opposite 
sign at the monolayers, thus indicating the different orientations of the molecules. 



6 - 



Q_ 

„ 4 





1 i 1 i 

\ 1 












l 
I 
l 


II 
I 
I 
I 




. . . 1 























119 



119.5 



120 

R 



120.5 



121 



FIG. 1. Density profiles for a spherical vesicle at equilibrium, with bulk densities p£ = 9 1CF 7 , and p$ = 3.4 1CF 14 . Dashed 
line: density of the majority component p A (R). Dot-dashed line: density of the minority component p B (R). Solid line: full 
density profile p(R). The interaction potential parameters are K = 10 and e = 1.4. The hard sphere diameter is used to provide 
the units of distance and density. 



As first step we study the phase diagram for planar membranes, taking the vesicle radius, R v , to infinity. In 
Figure 2 we present the line of stability for e = 1.4, with K being a function of the total density of surfactant in 
the bulk solution, p = p^ + p$ , for different values of the difference in the chemical potential of the two species, 
(3Ap — \og(pQ / Pq). In the region below each line, the membranes have Aft > 0, i.e. they are unstable and would 
dissolve in the water. In the upper region Af2 < 0, so that the planar membranes grow taking particles from the 
bath; then the bulk density is depleted to the coexistence value. The main qualitative change, with respect to the 
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pure surfactant membranes comes from the possibility that membranes have a spontaneous symmetry breaking, which 
gives them asymmetric composition </> ^ 0. This effect happens as a second-order phase transition, at positions on 
the coexistence curves (the places marked by the arrows in Figure 2) which help in stabilizing the membranes and 
lower down the lines in the figure. It is important to remark that this segregation is only possible when there is 
a reasonable amount of the two species, A and B, in the membrane. But the effect is obtained even if the bulk 
solution has strong majority of one of the species. The bilayer membrane is able to concentrate the minority species 
by many orders of magnitude of its bulk concentration, and then it may segregate preferentially in one of the the two 
molecular layers. Obviously, once the planar membrane undergoes the phase transition, there is no symmetry reason 
to avoid spontaneous curvature of the membrane into a closed vesicle. This is the coupling, between the composition 
asymmetry and spontaneous curvature, which acts to stabilize the vesicles proposed by Safran et al |10| , |l]| , and we 
may now analyze it within our microscopic model. The main difference here is that in our model this coupling is 
fixed by the molecular interactions and by the requirements of global thermodynamic stability, instead of being a free 
parameter (assumed constant) in the effective Landau free energy. 




log(p ) 

FIG. 2. Equilibrium lines in a log(po) _ K phase diagram for planar membranes, with e = 1.4. The dashed line gives the 
monocomponent limit A/i — + oo. The full lines show the coexistence line for several values of Afi — 4, 6, 8, from bottom to top. 
The arrows mark the points where <f> becomes different from zero. 

The first result of the consistency within our microscopic model, in which we have assumed perfect symmetry be- 
tween the two types of surfactant, is that the spontaneous curvature should vanish not only for symmetric membranes, 
with <j) = 0, but also for membranes with preferential segregation of each surfactant in one of the layers, 0^0, but 
with equal global concentration of the two species, ip = 0. Thus, in the Landau free energy the coupling term has 
to go like C(f>f(ip) (where f is an odd function of tp that vanishes at tp = 0), instead of being a constant times Ctfi, as 
taken in ref 11. In real catanionic mixtures, this symmetry between the two components is not perfect, but still they 
are probably closer to our ideally symmetric model than to the the model with constant f{ip)- 

The simplest Landau free energy we may write, with this requirement, is 

/ = Ti0 2 + T 2 ^ 2 + a l( j) A + a 2 <if> 4 + a 12 4> 2 ^ 2 + k[c 2 - 2 7 C<^] - 2^Afi, 

where we have approximated the function / by a linear function of ip. The coefficients a±, a 2 , a\ 2 , k and 7 are assumed 
to be positive constants. T\ and T 2 correspond to the reduced critical temperatures, which are negative below critical 
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points for the order parameters <f> and ip m the planar membranes (c = 0). In the effective free energy used by 
MacKintosh and Safran jy], the two temperatures Ti and T2 are taken to be identical, which would correspond to 
independent phase transitions in each monolayer, with the preferential segregation of each type of surfactants on 
uncorrelated patches in the two sides of the bilayer. In our case (with the term "fC(j)ip instead of "fCtp) we require 
T\ < T2, so that the parameter </> goes through a phase transition when T\ changes sign while Ti is always positive. In 
the case of catanionic surfactants our assumption is apparently easier to justify, given long ranged attraction between 
molecules of different types, than the hypothesis of MacKintosh and Safran, which requires an effective attraction 
between molecules of the same species in order to have segregation at each monolayer. 
The minimization with respect to the curvature leads to an effective free energy 

/ = T 1( f> 2 + T 2 i/j 2 + fll 4 + a 2 V 4 + (012 + 3k 7 2 )^V - 2(A M )V>. (6) 

There is a line of critical temperatures T c (A/x), starting at T c (0) = Ti, below which both C and cj> become nonzero. 
In Figure 3, we present a sketch of equilibrium curvature of vesicles, at fixed T, as a function of the Au in the free 
energy given by eq 6 — dot-dashed line — and by the free energy analyzed by MacKintosh and Safran — dashed 
line — . The qualitative difference is the behavior at Afi — associated to the symmetry of our model. In both cases, 
the general prediction of the Landau theory is that there may exist a tricritical point at which the transition becomes 
first order, through the coupling with ip. However, in our microscopic model, we have not found any situation in 
which the effective parameters in eq 6 lead to a first-order transition. The requirement of thermodynamic equilibrium 
with the bulk imply strong limitations in the possible values of these parameters. 




FIG. 3. Curvature of the vesicle C as a function the difference between chemical potentials Afi. Solid line: our microscopic 
theory along of the equilibrium curve, with K = 10 and e — 1.4. Dashed line: schematic curve obtained with a Landau theory 
assuming a coupling C(j> between the curvature and the membrane composition. Dot-dashed line: As for the previous curve, 
but now assuming a coupling Ccfit/j. 

If we compare the empirical form of the free energy given by the eq 6 with the numerical result for our microscopic 
model, we find that the fit of the free energy, as a function of curvature, to a second-order polynomial is always very 
accurate. Thus confirming a Helfrich-like picture in which an expansion of the free energy to two second-order in 
the curvature is carried out. The fit also allows us to safely extrapolate to the infinite radius (zero curvature) limit. 
However, the dependence of the free energy on the composition of the membranes, i.e. with the variables ip and <p, is 
quite beyond the truncated Taylor expansion used in eq 6. The actual results for the curvature, also shown in Figure 
3 — solid line — , show a richer dependence of the spontaneous curvature with A/i. 

To obtain the global phase diagram for curved structures in our model, we have taken as a starting point the results 
for the planar membrane C = 0, presented in Figure 2. In this way, for e = 1.4, we choose a value of the temperature, 
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for instance K = 10, so that for any value of A/i there is a unique value of the total bulk density of the surfactants, p , 
such that the planar membrane has zero surface tension, AQ(po, C = 0) = 0. When the dependence of O with respect 
to the vesicle curvature, C, is included we obtain results like those in Figure 4 for Ap = 18.718 (curve (a)) and 17.692 
(curve (b)). In the first case, the global minimum is precisely at c = and the membrane does not have spontaneous 
curvature. In case (b) the spontaneous curvature lowers the value of Ail below zero. The global thermodynamic 
equilibrium in this latter case requires a change of the bulk density po, from 1.1761 10~ 6 to 1.1756 10 -6 to recover the 
condition rnin[AQ(p, c)] = 0, as shown in curve (c). Now, repeating the procedure for any values of K, we could draw 
the phase diagram for curved structures (vesicles), but it is not necessary because is equal to the obtained for planar 
membranes since the shift of the bulk coexistence density is too small to be distinguishable in the scale of Figure 2. 
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FIG. 4. Excess free energy per particle, Au> — AQ/AN, versus curvature for K — 10 and e = 1.4. a) Planar membrane 
C = at coexistence, with log(po) = —13.5837 and Ap — 18.718; b) Vesicle C / beyond coexistence, with log(po) = —13.6533 
and A/i = 17.692; c) Vesicle C ^ at coexistence, with log(po) = —13.6537 and Ap — 17.692. Notice how slight the correction 
is between the parameters in b) and c) . 

In the Figure 5 we present another view of the coexistence line for bilayer vesicles, in terms of the concentration 
ratio and of the total bulk density (always for K = 10 and e = 1.4). The full line gives the concentration ratio 
measured in the bulk log(pQ / Pq) (which is precisely /3Ap), moving along a horizontal line would correspond to 
change the total concentration of surfactant but keeping the bulk ratio between the two species. Bulk dilutions on 
the left of the full line would be stable, while on the right of the line they would form bilayer membranes. This 
coexistence line also can be represented in term of the composition of the membrane instead of the bulk composition. 
So, in the same figure, the dashed lines give the ratio between the two surfactant species in the membranes, within 
our parametrization p^ = rjix\ + r\ix%. The ratio between the two species is directly given by the parameter ip, 
Pm/Pm — (1 — 2-0)/ (1 + 2%j)). As already noted for the planar case, the condensation of the surfactant on the 
membrane produces a very important change in the relative concentrations of the two surfactant species, given by the 
vertical separation between the full and the dashed lines in the Figure 5. A bulk dilution in which the species B is 
in clear minority (e.g. Ap = log(po / p^) ~ 10) coexists with vesicles in which the two components have comparable 
concentrations (log(p^/p^) ~ 1). 
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log(po) 

FIG. 5. Equilibrium curve for the vesicles as a function of the relative concentration and the total bulk density log(po). For 
the solid line, the relative concentration is measured in the bulk log(p^/p^) = Ap. Alternately for the dashed line, the relative 
concentration is measured in the membrane p^/Pm = (1/2 — tp)/(l/2 + tp). The energy parameters are K = 10 and e = 1.4. 
Notice how slowly the membrane ratio grows as compared with the bulk one, reflecting the ability of vesicles to concentrate 
the minority component. As explained in the text, the circles mark the places where C reaches zero. 

The results for the spontaneous curvature of the vesicles along the coexistence line are given by the full line in 
Figure 3, to be compared with the results of the simple Landau theories. For |Ap| > 18.7 the membrane is symmetric 
(and planar), since the minority component is too scarce to allow the segregation phase transition. This situation 
corresponds to the behavior on the right hand side of Figure 5, in which both the bulk and the membrane tend 
asymptotically to the one-component case. As we lower the values of |A/i| the transition takes place, there appears 
spontaneous curvature of the vesicles, and the composition asymmetry of the two layers develops rapidly. But this 
transition requires the incorporation into the membrane of a significant amount of the minority component, which 
pushes the membranes towards more symmetric composition, p^ « p^, and it is represented by the dashed lines going 
to zero, on the left of Figure 5. In accordance with the symmetry arguments presented above, the membranes become 
flat again at the symmetric mixture (Ap = and hence tp = 0). The three points marked by black circles along the 
coexistence line, in Figure 5, correspond to these transitions from planar to spherical vesicles. The existence of these 
transitions and the general shape of the phase diagram in Figure 5 is generic for symmetric surfactant mixtures and 
can be obtained from the simple Landau free energy (eq 6). What is more amazing in our microscopic calculation 
is the existence of intermediate points, at A/i « ±12, where the curvature also vanishes, although both cp and tp are 
nonzero. These points, marked by open circles in Figure 5, have unbalanced total composition both at the bulk and 
the membranes (po / Po ~ 10 5 , Pml Pm ~ 1-5) and the two layers of the membrane have broken the symmetry between 
the two species, but there still is no spontaneous curvature. This behavior is not predicted by simple Landau theories 
and implies that the coupling between <p and C has a complex dependence on tp. 

In order to allow better visualization of the composition of the bilayer vesicle, we show in Figure 6 the relative 
compositions of the A and B species for each of the layers. For high values of Ap > 18.7, in the flat membrane regime, 
we see the symmetry between the two layers, x\ = x\. As we decrease Ap below the transition to spherical vesicles, 
the curvature grows and the minority component B goes preferentially to the outer layer of the membrane (labeled 
as 1), but even in that layer the species B is still in minority. The curvature reaches a maximum for Ap; w 16, with 
vesicles of radius about one hundred molecular diameters. For still lower values of Ap the vesicles become larger 
until for Ap w 12 the spontaneous curvature disappears again. At this point the composition of the outer layer in 
the vesicles has become nearly symmetric, x^ w xf , while the inner layer keeps a clear majority of the species A. 
When we go to further lower values of Ap the membrane bends over in the opposite way, leaving now the layer with 
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majority of component B as the inner layer (indicated by negative curvatures in Figure 3 ). At A/i w 6.3 the size 
of the spontaneous vesicles reaches the minimum, with a radius of about twenty molecular diameters. The curvature 
goes again to zero if we approach the symmetric mixture Ap = 0, while the composition of the two layers is nearly 
frozen, with nearly complete segregation of each species to different layers. 




FIG. 6. The relative compositions x\ of each species in each layer versus A/i, along the coexistence line, for K = 10 and 
€ = 1.4. 



We have not found a direct explanation of this behavior, which may depend on the microscopic interactions of the 
model, but the general feature seems to be quite robust. The results for e = 1.2 (i.e. decreasing the difference between 
the A-A, B-B and the A-B interactions) show that the range of A/i with spontaneous curvature decreases but there 
is still a change of sign in C, with the same qualitative trend as for the case with e = 1.4. 

Making use of the information presented in Figure 5 it is possible to construct a phase diagram in terms of the 
usual experimental variables, the total densities of each amphiphile in the dissolution: p\ = p| + p^ n a m , where a m 
is the total amount of membrane area by unit volume of the dissolution. The two densities p^ and pf determine a 
point in the diagram; if this point is located on the left of the coexistence (solid) line, the bulk densities are equal to 
the total ones; that is, no structures are formed and a m = 0. If, on the other hand, the point lies on the right, the 
molecules can either stay on the bulk or form aggregates; the solution shifts towards a point on the coexistence line. 
The final equilibrium point can be easily determined from the knowledge of the composition of the aggregates, i.e. 
the dashed line. With this information we redraw the phase diagram, Figure 7, in the plane log(pt) - log(pf/pf ). 
We can distinguish three zones: one in which no aggregates are formed (I), another one in which they form planar 
membranes (II), and a third one where they form vesicles. The latter is divided in two subregions, corresponding to 
positive (Ilia) and negative (Illb) curvature. We believe these features of our model to be in qualitative agreement 
with experiments with catanionic systems in the dilute region ]l5|,p|,^6|~^8[ , specially the existence of two regions where 
vesicles are stable (zones HI) when the mixture is not equivalent, separated by a region in which the two components 
are present in similar amounts with flat membranes forming and precipitating. 
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log(p,) 

FIG. 7. Phase diagram in the log(pt) - log(p^/p^) (which is not Ap, since now we refer to total densities) plane, at fixed 
K = 10 and e = 1.4. The solid line is the coexistence line between free molecules and aggregates; the dashed ones, between 
different types of aggregates. The meaning of the regions and the dot-dashed paths is explained in the text. The diagram is 
completely up-down symmetric, but for the sake of clarity we only draw the regions in the upper side and the different paths 
in the lower one. 



In order to compare our results with experiments which measure the evolution of the average radius of the vesicle 
when water is added to the solution 0, we have consider the behavior of our system when the ratio of the two 
components, pf/pf, is kept fixed while the total density p t — pf + pf is varied. This process would follow horizontal 
lines in Figure 7. Depending of the value of \og{pf / pf) we can have different behavior; we show three interesting 
examples, labeled from 1 to 3. The sequence of structures, from low to high values of p t , would be: for the first 
case, no structures-planar membranes; for the second, no structures-vesicles-planar membranes; and for the third, 
no structures-vesicles. 
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FIG. 8. Average vesicle radius as a function of the total density, keeping the composition ratio of the two components fixed. 
In panel (a) we present the results obtained with our model, following two horizontal paths of type 3 in Figure 7, as a function 
of log(pt). The full line corresponds to log(pf /p^) = —1, and the dashed line to —0.1. Panel (b) show the experimental results 
of Kaler et al in ref 2, for composition ratios 1/9 (solid line) and 1/4 (dashed line). The radii are given in nanometers and the 
total density is measured by log(l — w/100), where w is the water-weight percentage. 



The evolution of the radius with pt depends strongly on the ratio / pf and it may be nonmonotonic. In Figure 8(a) 
we present our results for two different cases, both corresponding to a sequence of structure type 3. The nonmonotonic 
dependence comes directly from the nonmonotonic behavior of C in Figure 3; it would not appear in a Landau free 
energy analysis with a constant coupling between C and <j). The experimental observation of the spontaneous curvature 
presents strong technical difficulties, with uncertainties in the equilibration and the polydispersity. However, the 
results, for mixtures of cetyl trimethilammonium tosylate and sodium dodecylbenzene sulfonate ||, shown in the panel 
(b) of the Figure 8, present clear examples of non-monotonic dependence of the radius with the total concentration 
of surfactant, p t . The fragmentary experimental information does not allow us to confirm or to refute the presence of 
the intermediate point with C = found in our model. A direct comparison with our predictions would require the 
chemical analysis of the inner and the outer layers of the vesicles, which has not been carried out so far. As for the 
typical size of the spontaneous vesicles, we can only expect a qualitative agreement with our simplified model. In our 
results, the vesicles have always radius larger than 20 and typically in the range of 100 in units of the molecular size 
(the hard core parameter); the experimental values ranging between 20 and 80 nm are in the same order of magnitude, 
for molecular sizes in the nanometer range. 



IV. CONCLUSIONS 



We have developed a self-consistent microscopic model for the study of aggregates formed in a mixture of two 
different amphiphilic species (A and B) in water. The simplicity of the interactions considered here makes possible a 
self-consistent thermodynamic treatment, assuring thermodynamic equilibrium between the solution of amphiphilcs 
and the aggregates. We have assumed that the interactions between the amphiphiles are symmetric, the A— A in- 
teraction is equal to the B-B one but different from the A-B interaction. Using a semiempirical Landau theory, 
we — as Safran et al did previously p0| , pl| — have found that the mixture of interacting amphiphiles stabilizes the 
spherical vesicle with respect to planar membranes. However, there are important differences between the two models: 
Safran et al use a Landau free energy based on the hypothesis that the driving force for the spontaneous curvature 
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is the existence of independent phase transitions in each monolayer, due to an effective intralayer repulsion between 
the molecules of different types. In our case, we assume that the driving force is an interlayer attraction between 
different molecules in opposite layers of the membrane. Our hypothesis seems to be easier to justify for real catanionic 
mixtures and leads to a different coupling between the asymmetry of the molecular composition and the curvature of 
the vesicles. In our case, the coupling term go like C<pf(ip), where f is an odd function, thus it vanishes not only for 
membranes with equal composition in the two layers ((/> = 0) but also for membranes with different layers <j> ^ but 
with equal global compositions of the two surfactants (tp = 0). 

Moreover, in our density functional approach, we obtain the free energy directly from molecular interactions, 
instead of having empirical parameters as in any Landau free energy. Our results show that the presence of first- 
order phase transitions, predicted by the different Landau theories, may be easily frustrated by the requirements of 
global thermodynamic stability. The membranes may be regarded as two dimensional phases, but they have to be at 
coexistence with a diluted bulk and, in our microscopic model, this imposes severe restrictions in the possible behavior 
of the free energy. In fact, the truncated expansion of the free energy in terms of the order parameters <p and tp turns 
out not no be so useful, since the coefficients vary in a quite complicated manner as the thermodynamic conditions 
are varied — despite the truncation being, of course, correct in a limited range of application, i.e. small values of the 
order parameters. 

The phase diagram obtained with our model is in qualitative agreement with the experimental diagram for a 
system comprised of an anionic and a cationic surfactant in the dilute region. We have also studied the evolution of 
the vesicles size when water is added to the dissolution. We have found that the radius of the vesicles increase or 
decrease depending of the ratio p A /p B . This fact could explain the complicated experimental evolution of the vesicle 
radius, although experimental technical difficulties can also have an influence. 

In conclusion, the minimal model for the mixed anionic-cationic surfactant vesicles, presented in this work, is a 
fundamental description of the thermodynamic of these systems. Although the micellar phase has not been taken into 
account in the present work, this fact does not change the main conclusions since the vesicles and plane membranes 
are always more stable than the micelles. This assumption is based on the fact that we have used q = 0.5 for 
all the intermolecular potentials, and this value stabilizes the plane membranes compared to the micelles for the 
monocomponent case. Also, we have neglected the entropy of shape fluctuations and the center of mass translations 
for the aggregates. The inclusion of such terms for vesicles would give a finite width to the 'coexistence line' calculated 
here. However, given the large size of the spontaneous vesicles (always with more than 10 4 molecules) and the large 
value of the bending rigidity, the effects would be negligible. Extension of the present model to the study of mixtures 
of amphiphiles with different self-assembly behavior, including micellar aggregates, is being in progress. 
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